UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



LIBRARY ROUTINE M ^h - 318 



TITLE: 

TYPE; 

CAPACITY: 

ACCURACY: 

DESCRIPTION: 



METHOD OF USE: 



Symmetric Matrix Inversion 

Entire program 

2 < n < 139 where n is the order of the matrix 

A function of the order and conditioning of the matrix 

This routine inputs a symmetric matrix as a lower triangular 

matrix, and inverts it. Input may be from tape or from the 

drum. Output may be stored on the drum as a lower triangular 

matrix j or punched to "d" decimal places (l < d < 12) as a 

square matrix, a lower triangular matrix, or a column of the 

diagonal entries. 

STOPS 

Master tape 3022K 

Parameter tape 2^068 

Data tape 2022K 

At the end of a problem the program will stop on a 2022K order, 
Raising the black switch will begin a new problem with new 
parameters • Moving the white switch up and down will begin 
a new problem with the old parameters* 
THE PARAMETER AND DATA TAPES: 

The parameter tape is punched as+n+s+d+pT 

i) n is the order of the matrix 

ii) s is some integer power of 10, This is used "to scale the 

input matrix so that the largest element in the matrix is as 

close to l/2 in machine representation as possible, without 

equalling or exceeding it. ; 

iii) d is the number of decimal places in the output 

iv) p is the input option 

a) designates tape input. It need not be punched. 

b) any non-zero number designates drum input beginning 
at location 2906. 
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EXAMPLE 1; 



v) T is the output option 

a) N designates punching a lower triangular matrix 

b) J designates punching a square matrix 

c) F designates punching a column of the diagonal matrix 

d) L designates 'storing q£..;the lower triangular matrix 
on the drum beginning at drum location 2906, with a 
sum check at the end of each row. 

The data tape is punched as signed fractions by rows . It is 
in the form of a lower triangular matrix and the matrix is 
terminated by a J. 

Input the following matrix from tape and punch the inverted 
matrix as a square matrix to 5 decimal places 



.5 -.6 -S 
-.6 .8 .9 



EXAMPLE 2: 



EXAMPLE 32 



The largest element would be greater than l/2 in machine re- 
presentation. Scale the matrix by 10*1, that is s = -1. The 
.parameter tape is punched a* +3 -1 -+5-1-+J. ] The data tape, is punched 
as +5 : ..--*6 .-+8l-.^6-:-h-9 -9J,-i- - '; «. 

The matrix in example 1 is already stored on the drum as a 
lower triangular matrix, beginning at location 2906* Invert 
and punch to 6 decimal places as a column of diagonal entries • 
The parameter tape is punched as +3 -1 +6 -kLF. 
Invert the following matrix and store on the drum* 



T 








,2 ^0J 





v)0 .3 





.k 



The largest element is less than 1/2 in machine representation. 
Thus s = 0, Notes also that d = since there is no tape out- 
put „ The parameter tape is punched as +4L. The data tape is 
punched as +1 + +2 ++ +3 +++ +4 J- 
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THE OUTPUT MATRIX: 

The elements of the punched output matrix are punched as signed 

fractions , There are four modes of output • 

i) Square matrix . The elements are punched by rows with an 

N after each row and a J after the matrix, 
ii) Lower triangular matrix * The elements are punched by- 
rows and 2 line feeds and carriage returns follow each row* 
A J will follow the matrix * 
iii) Column of diagonal entries * Eachcelement' is followed by 

a line feed and carriage return, 
iv) The output is stored as a lower triangular matrix on the 
drum beginning at location 2906. A sum check is stored 
after each row. 
The punched outputs have a d-place number appearing after the 
output matrix. This indicates the position of the decimal 
point as lying after the column in which a 1 appears. 
EXAMPLE k: The output to example 1 appears as 

+18889 +13333 +007*aN 
+13333 +10000 +01111N 
+007+1 +01111 -OO^NJ 
+01000 
This would be interpreted as 

T8.889 13.333 0,7^1" 
13.333 10.000 J 1.111 
_0.7+l " .1.111 -0.^-_ 

In general: Let A be the matrix (of fractions) punched on the 
input tape and let B be the output matrix (with decimal points 
in positions as specified above). 
Thens 

AB = I(10" S_1 ) 
where I is the unit matrix and s is the input s-cale factor dis- 
cussed above. Hence 

A" 1 = 10 S+1 B 
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DURATION: It takes a minute to read in the master tape. The computation 

time is given "by fb.OO^n + 0.08n + O.lnJ seconds. 
PROGRAMMED STOPS ; 

FFOOO Matrix is singular 

FFOOl Matrix is overs caled 

FF002 Diagonal element too small 

FF003 Master tape read incorrectly. Reread. 

Moving the white switch up and down will start a new problem 

at stops FFOOO to FF002, 
MATHEMATICAL METHOD 1 

Let A "be the matrix to be inverted. Let L be a lower, U an 

upper triangular matrix, and D a diagonal matrix. Then A = LDU, 

If A is symmetric A = (LDU) T = U T P T L T » 

That is L = U T , D = D T and U - L T 

A" 1 = (LDU)" 1 s U _1 D" 1 L _1 = (l t ) _1 d" 1 L" 1 . (*) 

The program stores A as a lower triangular matrix on the drum. 

By a series of row transformation on A, L, and D are obtained. 

L~ is obtained by the following; 

if j < k 



V. 



-Z J . - l / . if j > k 

m=j + l ^ mk 



D" is obtained by the following: 







A is obtained from equation (*). 
A necessary condition for singular 
machine representation, due to the accumulation of round-off 



A necessary condition for singularity is that d. . = 0, In 

J J 
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errors, it is difficult to detect a near singular matrix. The 
results should always be checked when singularity is detected 
by the program, 



nj 



DATE March 20, 


196l 


PROGRAMMED BY F 


, Shimamoto \ 


APPROVED BY \)riAAv"\A2^ 





LOCATION 



10 
11 
12 
13 
14 
15 
16 

17 
18 

19 
20 
21 



ORDER 



(N) 

(S) 

(P) 

(I) 

(0) 

(MS) 

(Z) 

(VI) 

(V2) 

(V3) 

(1) 
(1-D 



003K 
00 F 
00 607F 
00 F 
00 7^6F 
00 F 
00 885F 
00 F 

00 2906F 

0010K 
00 F .> 
00 F 
00 F 
00 F 
00 F 
00 F 
00 F 
00 F 
00 F 
00 F 
00 F 
00 F 
00 F 
00 F 
7J S 
3^ 0S3 
00 F 
00 Qk 
Ik S5 
kO S5 
00 F 
00 IF 
00 IF 
00 IF 



7 



NOTES 



PAGE 1 M 3^ 



Location of vectors A, B, C 



Location of matrix on drum 



Constants 



LOCATION 


ORDER 


NOTES PAGE 2 


m y\ 


22 


(+) 


00 F 
00 10F 








23 


(-) 


00 F 












00 1000 0000 0000 J 






2k 


(+D 


00 F j 

00 9999 9999 9999 J 






25 


w 


00 F 
00 3^F 








26 


(ci) 


00 F 
00 F 








27 




00 F 
00 2F 








28 


(C2) 


26 (Yl) 












00 2S6 










005^5K 






5^5 


(mtr) 


JO 27F 
50 L 


Store part 1 of program on drum 




5^6 




26 (Yl) 

00 2560F 






5W 




00 h6f 
26 999F 






5^8 




JO 29F 
50 3L 






5^9 




26 (Yl) 

00 2677F 


Store part 2 




550 




00 125F 
26 999F 






551 




JO 52F 
50 6l 


4 




552 




26 (Yl) 

00 2803F 


► Store part 3 




553 




00 102F 
26 999F 


j 




55^ 




50 27F 


1 




"' 




50 9L 


1 





LOCATION 


ORDER 


NOTES PAGE 3 M 3^ 


555 


26 (Yl) 


f 




00 2560F 


J Playback part 1 


556 


00 116F 
22.29F 




557 


50 29F 






50 12L 


r Playback part 2 


558 


26 (.Yl) 

00 2677F 




559 


00 125F 






26 29F J 


560 


50 52F 


) 




50 15L 


k Playback part 3 


561 


26 (Yl) 

00 2803F 




562 


00 102F J 




26 53F 




563 


20 9L 






22 19L 


End of program; bl. sw. to read new para- 


56^ 


50 27F 
50 19L 


meter, wh. sw. to retain old parameter 


565 


26 (Yl) 


1 




00 2560F 


r Playback part 1 


566 


00 ii6f 
26 31F 


J 


". 56? 


(Yl), OOK 


Drum transfer routine (Yl) 

\ 




00 29K 


29 


(IN) 52 10F 


] 




50 L 


> Read in parameters 


30 


26 (N12) 

ko i4f 


I 


31 


L3 (S) 
36 l^L 




32 


L5 (S) 
32 8L 




33 


kO kF 
L5 M 





LOCATION 


ORDER 
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3^ 


40 3F 






F5 4f 


Store negative scaler 


35 


36 14L 
40 4F 




36 


50 3F 
7J (m) 




37 


26 5L 
LI (S) 




38 


4o 4f 
L5 (+) 




39 


40 3F 






F5 4f 


Store positive scaler 


40 


36 14L 
4o 4f 




41 


50 3F 
75 (+) 




42 


S3 F 
26 10L 




43 


L5 (1) 

40 5F 




44 


LI (N) 






40 4F 




45 


L3 (3E) 






32 46L 


Input option test: jump to tape input 


46 


L3 (S) 
32 64l 


1 


47 


50 S3 


N 




50 18L 


Playback row of matrix from drum 


48 


26 (Yl) 

00 s6 




49 


00 IF 




50 


50 20L 

26 (m) 

22 22L 


J 


51 


JO S3 
50 22L 


/ ' . ' - ■ ......... 



LOCATION 


ORDER 


NOTES PAGE 5 M 34 


52 


26 (Yl) 


Store row of matrix on drum 




00 S6 


53 


00 IF 
F5 4F 




5^ 


32 64L 
40 4F 




55 


F5 5F 
40 5F 




56 


L4 19L 
40 19L 




57 


40 23L 
L5 20L 




58 


L4 (1-1) 
46 20L 




59 


46 24L 
22 18L 




60 


$a) K5 f 


} 




42 38L 


1 +s 

1 Scale a row of matrix by 10— routine 


6l 


42 43L 
LI 5F 




62 


4o 6f 
LI (S) 




63 


32 ^OL 
L5 (VI) 




64 


42 36L 
42 37L 




65 


50 3F 
75 F 




66 


00 39F 
40 F 




67 


F5V6F 
32 F 




68 


4o 6f 
F5 36L 




69 


26 35L 






L5 (VI) 



LOCATION 


ORDER 
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70 


kO 42L 
50 3F 




71 


7J F 
kO F 




72 


F5 6f 
32 F 




73 


ko 6f 

L5 42L 




7 4 


14 (1-1) 






26 4lL 


<L 


75 


50 F 

4i 6f 




76 


L5 74l 






40 23(N12) 


\ 


77 


50 S3 


) 




50 48L 


[ Read In data 


78 


24 ($02) 






F5, $F\ 


; 


79 


4o 6f 
LO 5F 




80 


32 52L 
L^ 2F 




81 


.26 4(N12) 
L3 (S) 




82 


32 55L 
50 53L 




83 


26 (m) 


1 




22 55L 




84 


JO S3 


N 




50 55L 




85 


26 (Yl) 






00 S6 J / Store a row of A on drum 1 


. 86 


00 IF 1 
F5 kF I 




87 


32 64L I: 1 




40 4f I , 

1 , : ■■ \Mm , , , , { 



LOCATION 



ORDER 



NOTES 
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WW 



88 

89 

90 

91 
92 

93 
9 k 

95 
96 
97 
98 

99 

100 

101 

102 

103 

\0k 





F5'5E 
kO 5F 

14 56L 

40 56L 
L5 57L 
IA (l-l) 
46 57L 
L5 (VI) 

if6 2l(N12) 

hi 6f 

22 51L 

L5 (N) 

l4 66l 

^0 5F 

50 (1) 

41 S5 
F5 66l 

42 66l 
LO 5F 

32 69L 

22 66l 

L5 (N) 

00 20F 

l4 (N) 

kO 193F 

L5 (1) 

ko 256F 

40 2F 

41 (MS) 

22 12(MTR) 
22 49 (IN) 
00 F 
(N12 )00K 
OOK 
L3 F 
kO 607F 



Clear C vector 



Input routine (N12) 



LOCATION 


, 1 

ORDER 


NOTES PAGE 8 M 34 


1 


22 5^5F 
50 F 
26 L 
26 IN 

0029K 




29 


F5 (VI) 
42 257F 




30 


F5 (1) 
40 3F 




31 


LI (N) 
L4 (1) 




32 


40 j(Cl) 

4o 7F 




33 


L5 (l) 






4o 6f 




34 


85 11F 


■ 




00 s6 




35 


4o s4 






L5 S5 




36 


40 4f 
22. 8L 


; 


37 


50 S3 


-) 




50 8L 


! -. 

j w Playback one row of matrix 


38 


26 (Yl) 


1 




00 2S6 


39 


00 2F 






LI 1S5 




40 


L4 4f 






40 IF 




41 


L5 S3 






40 1S3 


! 


42 


L3 IF 


! 




32 25L 


i 


^3 


L7 IF 






00 20F { 


,,... „ . , _iJ 



LOCATION 



ORDER 



NOTES 
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44 


46 20L 




46 24l 


45 


Ll IF 




32 23L 


46 


L4 4f 




40 4f 


^7 


L5 (V2) 




42 19L 


48 


42 20L 




L5 F 


49 


10 F 




40 F 


50 


F5 20L 




LO 24L 


51 


36 25L 




F5 19L 


52 


22 18L 




L5 1S4 


53 


10 F 




40 1S4 


54 


41 IF 




L5 S3 


55 


40 F 




L7 S4 


56 


12 F 




32 30L 


57 


F5 IF 




40 IF 


58 


L5 F 




10 IF 


59 


26 26L 




L3 IF 


6o 


32 39L 




L7 13 1 


6i 


00 20F 




46 36L 



Test: Is a - a Jk < 0? 



LOCATION 


ORDER 


NOTES PAGE 10 M 34 


62 


Fl 2F 
4o 5F 




63 


L5 (VI) 

k2 35L 




64 


k2 36L 






L5 F 




65 


10 F 
kO F 




66 


F5 5F 
kO 5F 




67 


32 39L 
F5 35L 




68 


22 3^L 
L7 kF 




69 


l4 (1) 






00 20F 




70 


46 k^L 
L3 S4 




71 


36 124L 
50 S3 


■\ 


72 


75 (+1) 






66 sk 


r Form a ., 


73 


S5 F 
4o 8f 

i 


a. . ■ 
JO 

! 


74 


10 IF ' 

00 u 1 ! 


i 


75 


32 46l 


■ 




kO S3 


1 


76 


LI 6f 
'40 5F 




77 


F5 43L 
42 50L 


J 


78 


F5 42L 
42 51L 


i 

I 


79 


42 54L 


j 


i 


L5 F 


i 



LOCATION 


ORDER 


NOTES PAGE 11 M 5^ 


80 


kO F 
L5 F 




81 


kO 9F 
56 55L 




82 


50 8F 


1 




7$. ; F 


f Form /jk 


85 


IX 9F 
iK).F 


J 


8k 


F5 5F 
kO 5F 




85 


52 58L 
F5 50L 




86 


k-2 50L 
F5 51L 




87 


22 k9h 
F3i 2F 




88 


40 F 
L5 (VI) 




89 


^2 6lL 
42 62L 




90 


L7 5F 




91 


52 65L 
L7 F 




92 


40 5F 
F5 F 




95 


40 F 
56 66L 




9k 


F5 6lL 






26 60L 




95 


Fl 2F 






kO\ 9F 




96 


LL 5F 
52 "11 "IT, 




97 


F5 F 






\2 71k 





LOCATION 


ORDER 


NOTES PAGE 12 M 3k 


98 


FO (l) 

LO IF 




99 


kO F 

L3 F 




100 


32 jSl 

L5 (VI) 




101 


k2 73l 
h2 75L 




102 


50 (1) 

L5 F 




103 


10 2F 
00 IF 




iafr 


32 75L 
kO F 




105 


F5 9F 
32 78l" 




106 


kO 9F 
F5 73L 




107 


26 72L 
L5 is 5 




108 


lk F 
ko 1S5 




109 


JO S3 


\ 




50 80L 


L Record a row of L on drum 


110 


26 (Yl) 
00 2S6 




111 


00 2F 






F5 7F 


/ 


112 


32 93L 
ko 7F 




. 113 


L5 10L 

IA (1-1) 




llii- 


kO 10L 
k6 82L 




115 


42 78L 






^2 79L 





LOCATION 


ORDER 


116 


F5 2F 




40 2F 


117 


L4 (1) 




l4 9L 


118 


40 9L 




40 81L 


119 


F5 12L 




42 12L 


120 


42 23L 




42 24L 


121 


F5 6f 




4q 6f 


122 


22 8L 


j 


F5 (Cl) 


123* 


36 115L 




F5 6l 


124 


42 6l 




l4 (1) 


125 


42 10L 




42 78L 


126 


42 79L 




F5 <va) 


127 


k2 12L 




42 23L 


128 


^4-2 24L 




L5 3F 


129 


40 2F 




l4 (1) 


130 


40 3F 




L4 (C2) 


131 


40 (C2) 




40 9L 


132 


40 8lL 




L5 3? 


133 


00 20F 




46 10L 



NOTES 



PAGE 13 M 34 



LOCATION 


ORDER 


NOTES PAGE 14 M 34"~~ 


13^ 


46 82L 
F5 25L 




135 


42 25L 
42 \2L 




136 


4£ kSL 
L5 3F 




137 


l4 5l 

40 % 




138 


L5 12L 

L4 (1-1) 




139 


46 12L 
F5 (Cl) 




l40 


26 3L 
F5 F 




i4i 


40 F 
L7 5F 




l42 


50 F 
00 IF 




3A3 


40 5F 
26 67L 




l44 


L5 (1) 
40 3F 




l45 


LI (N) 

40 F 




146 


40 5F 
L5 (V3) 




1A7 


42 119L 
42 120L 




l48 


L5 (MS) 

L4 F 




149 


32 121L 
L7 F 




150 


40 (MS) 

F5 F 




151 


32 15(MER) 






40 F 


- 



LOCATION 


j ORDER 


NOTES PAGE 15 M 34 


152 


F5 11% 
26 118L 




153 


FF F 






26 18(MER) 
OOK 


FFOOO; matrix is singular, vfcu sw, to 
begin new problem 





L3 F 
40 608F 




1 


22 548F 
50 F 
26 L 

26 in 

0052K 




52 


00 F 
00 66F 


>| 


53 


L5 196F 






40 3^9F 


Initialize addresses for rest of program 


5k 


L5 (VI) 




42 I76F 




55 


42 194F 
42 195F 




% 


42 268F 
42 279F 




57 


L5 (V2) 
U2 210F 




58 


\2 211F 
42 262F 




59 


42 264F 
42 280F 




60 


,42 281F 




- 


L5 8lF 




61 


40 246F 
40 410F 




62 


40 527F 
40 238F 




63 


L5 553F 






L4 196F 





LOCATION 


ORDER 


NOTES PAGE 16 M 34 


6k 


40 399F 
L5 (l-l) 




65 


46 312F 
k6 400F 




66 


46 4llF 
46 239? 




61 


-46 247F 
L5 (V3) 




68 


42 200F 
42 201F 




69 


42 233F 
40 159F 




70 


L5 138F 
4o 521F 




71 


40 546f 
4o 503F 




72 


F5 (M&) 
00 20F 

0073K 




73 


46 1L 






L5 (+1) 


Res c ale L routine 


74 


10 F 
00 IF 




75 


4o 6f 

L5 (VI) 




16 


42 12L 
42 13L 




11 


L5 3F 
L4 12L 




78 


40 4F 
L5 S5 




79 


L4 (MS) 
40 2F 




80 


50 S3 


) 




50 7L 


I Playback a row of L 



LOCATION 


ORDER 


NOTES PAGE 17 M 3k 


81 


26 (Yl) 
00 s6 


J 


82 


00 IF 
L3 2F 




83 


36 17L 

L7 2F 




& 


00 20F 
k6 13L 




85 


01 IF 


1 




L5 F 


V Rescale a row of L 


86 


10 F 
1*0 F 


1 


8T 


F5 12L 
k2 12L 




88 


hZ13L 
LO k¥ 




89 


36 17L 
22 12L 




90 


L5 S3 
ko S^ 




91 


L5 6f 
ko S3 


^% 


92 


JO S3 


] 




50 19L 


r Record a row of L 


93 


26 (Yl) 
00 s6 




9 k 


00 IF 






F5 5F 


J 


95 


32 31L 
kO 5F 




96 


F5 5L 
k2 5L 




97 


F5 3F 

^0 3F 




- 98 


lA 8l 






ko 8l 





LOCATION 
99 

100 

101 

102 

103 

10k 

105 

106 
107 
108 
109 

110 

111 

112 

113 
Ilk 

115 



(IT) 



ORDER 


NOTES PAGE 18 


M 3^ 


kO 20L 






L5 9L 






Lk '(1-1) 






k6 9L. 






k6 21X 






L5 1TL 






L4 (1-1) 






40 17L 






F5 i8l 






k2 18L 






22 2L 






M (z) 






LI (M3) 






40 S5 






OOK 






L5 (C2) 






ko 18L 






to 77L 






LI (I) 






il-0 2F 






I* (1) 






^0 IF 






L5 87L 






LO (1-1) 






ko 3F 






LO (1-1) 






Lk 22L 






^Q 22L 






LI IF 






Lk 88L 






^0 88L 






L5 70L 






LO 2F 






ko 87L 






L5 3F 
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NOTES PAGE 19 M 3^ 


116 


1* 53L 
ko 53L 




117 


k2 16L 
LI 2F 




118 


ko 3F 
kO ^F 




119 


00 20F 
k6 1% 




120 


k6 78L 






F5 (MS) 




121 


00 20F 
k6 30L 




122 


LI (MB) 
kO F 




123 


50 S3 


1 




50 17L 


" Playback row of L 


12^ 


26 (Yl) 
00 F 




125 


00 F 






L5 18L 


I 


126 


LO 91L 
L^ 3$kF 




127 


kO 5F 

ko 6f 




128 


LI S3 

^0 S3 




129 


F5 IF 
32 7<?L 




130 


L5 22L 
LO (1) 




131 


k2 29L 
lj-2 3l*-L 




132 


42 57L 
L.5 86l 




133 


^0 37L 






kl F 


< i " '■ ■ .111 ■ . 



LOCATION 


ORDER 


NOTES PAGE 20 M 34 


13^ 


4l 7F 
41 8F 




135 


50 7F 

L5 F 




136 


10 F 
00 IF 




137 


26 35L 
kO 3F 


■*v 


138 


85 11F 
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